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Abstract 



The lifetimes of the B + and B° mesons, and their ratio, have been measured in the OPAL 
experiment using 2.4 million hadronic Z° decays recorded at LEP. Z° — > bb decays were tagged 
ON ' using displaced secondary vertices and high momentum electrons and muons. The lifetimes were 

on : then measured using well-reconstructed charged and neutral secondary vertices selected in this 

X| tagged data sample. The results are 



r B+ = 1.643 ± 0.037 ± 0.025 ps 



D. r B o = 1.523 ± 0.057 ± 0.053 ps 

r B +/r B o = 1.079 ±0.064 ±0.041, 



where in each case the first error is statistical and the second systematic. 

A larger data sample of 3.1 million hadronic Z° decays has been used to search for CP and CPT 
violating effects by comparison of inclusive b and b hadron decays. No evidence for such effects is 
seen. The CP violation parameter Re(es) is measured to be 

Re(e B ) = 0.001 ± 0.014 ± 0.003 

and the fractional difference between b and b hadron lifetimes is measured to be 

Ar\ _ r(b hadron) — r(b hadron) 



r 



t( aver age) 



-0.001 ±0.012 ±0.008. 
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1 Introduction 



The lifetimes of b hadrons depend both on the strength of the b quark coupling to the lighter c and 
u quarks, and on the dynamics of b hadron decay. In the spectator model of heavy hadron decay, the 
decay of the heavy quark is unaffected by the presence of the other light quarks in the hadron, so the 
lifetimes of all hadrons containing the same heavy quark are predicted to be equal. This model fails 
badly for the charm hadrons, where the lifetime of the D + is more than 2.5 times the lifetime of the 
D° Q, but is a better approximation for the b hadrons, due to the larger mass of the b quark ||]. 
The difference between the B + and B° lifetimes^] is expected to be at most 10 %, and depends on the 
details of the various non-spectator processes contributing to their decay. Measurements at the level 
of a few percent or better are therefore needed to test this prediction and probe the dynamics of b 
hadron decays. 

Experimentally, most measurements of the B + and B° lifetimes have been performed using semilep- 
tonic decays, fully or partially reconstructing the decay products to distinguish B + from B° [||,|| . These 
measurements are limited due to the small branching ratios and limited reconstruction efficiencies for 
the selected final states. A more inclusive approach is to reconstruct resolvable secondary vertices from 
b hadron decays, since their long lifetimes lead to significant decay lengths In these analyses, 

the B + and B° decays are distinguished by reconstructing the charge of the secondary vertex. This 
technique results in much larger data samples, and offers the best hope of improving the precision. A 
measurement of this type is presented in the first part of this paper. 

Inclusive samples of B° decays can also be used to search for CP and CPT violating effects. 
Although CP violation has so far been seen only in the neutral kaon system, possibly large effects are 
expected also in the B-meson system, so it is worthwhile to search for them even with the relatively 
small data samples collected at LEP. In the B°-B° system the weak eigenstates |B°) and [B°) differ 
from the mass eigenstates |Bi) and IB2): 

| R v (l + 6 B +^)|B ) + (l-6B-^)iB°) 



|B 2 



(l_+€g -5 B )\B°) -(l-e B + 5 B )\B 
^/2(l + \e B -5 B [ 2 ) 



where the parameters e B and 5b parameterise indirect CP violation and CPT violation respectively J?j]. 
These parameters can be studied using semi-leptonic b hadron decays, and limits of a few 10 -2 have 



been set on both quantities ||q-|ll|. A non-zero value of e B is also expected to produce time dependent 
asymmetries in inclusive B° vs. inclusive B° decays . This provides a second method to look for CP 
violating effects in b decays, and such an asymmetry is searched for in the second part of this paper. 
The same data sample is also used to test a basic prediction of CPT invariance, that the lifetimes of 
b and b hadrons are equal. 

A brief overview of the analysis strategies is presented in Section |2[ followed by a review of the 
data and Monte Carlo samples in Section ^. The parts common to both B + and B° lifetime and 
CP(T) violation analyses, namely the bb event tagging and b hadron production flavour tagging, are 
discussed in Sections |] and [|. The lifetime analysis is described in detail in Section ^ and the CP 
and CPT violation analyses in Section |^. Finally, all the results and conclusions are summarised in 
Section @. 



2 Analysis Overview 

The analyses exploit the characteristic topology of the Z° — > bb decay: two back-to-back jets aligned 
along the direction of the thrust axis. The event is divided into two hemispheres by the plane per- 
pendicular to the thrust axis and containing the e + e~ interaction point. One hemisphere (the 'tag 

1 Charge conjugate states are implied when discussing the lifetimes of individual b hadron species. 
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hemisphere' T) is tagged as containing a b decay using either a displaced secondary vertex or a 
high momentum lepton. The production flavour of the b hadron in the tag hemisphere (whether it 
originated from a b quark or b antiquark) is also determined, using jet, vertex and lepton charge 
information. The unbiased b decay in the other hemisphere (the 'measurement hemisphere' M) is 
used to perform the measurement of b hadron decay time. 

The decay time of the b hadron in the measurement hemisphere is determined by reconstructing 
its decay vertex and energy. In the B + and B° lifetime analysis, this decay vertex is required to 
be significantly separated from the event primary vertex, and very strict requirements are placed on 
the confidence with which tracks are assigned to either this secondary vertex or the primary vertex. 
These requirements lead to a final data sample of only about 10 000 reconstructed vertices with well 
determined charge. The reconstructed decay times of the b hadrons giving rise to these charged and 
neutral vertices are then used to determine the lifetimes of the B + and B° mesons, employing the 
excess decay length technique []|, |L|] to eliminate biases caused by the separated vertex requirement. 
The correlation of the sign of reconstructed charged b hadrons in the measurement hemisphere with 
the production flavour of the b hadron in the tag hemisphere is used to measure the probability of 
mis-reconstructing the b hadron charge. 

In the CP(T) violation analysis, only the decay time of the b hadron in the measurement hemi- 
sphere is reconstructed. The charge is not determined, so no strict requirements are placed on the 
quality or separation of the b hadron vertex, leading to a much larger data sample of about 400 000 
events. The production flavour of this b hadron is inferred from that of the b hadron in the tag 
hemisphere, aided by information in the measurement hemisphere, and taking into account the effects 
of B° and B s mixing. The decay time distribution of decays tagged as b or b hadrons is then used to 
search for CP and CPT violating effects. 



3 Data sample and event simulation 



The OPAL detector is well described elsewhere [14-16]. The analyses described here rely mainly on 
charged particle track reconstruction using the central tracking chambers and the silicon microvertex 
detector. The b hadron lifetime analysis requires excellent pattern recognition, and uses only data 
taken between 1993 and 1995 with the upgraded silicon microvertex detector measuring tracks in both 
the r-cj) and r-z planes^ |l6fl . The CP(T) violation analysis also uses data taken in 1991 and 1992 with 
the original silicon microvertex detector measuring only in the r-(j) plane [15|. 

Hadronic Z° decays were selected using standard criteria [17], additionally requiring at least 



7 charged tracks in the event. The thrust axis polar angle cos#t was calculated using charged tracks 
and electromagnetic clusters not associated to any track. To ensure the event was well contained 
within the acceptance of the detector, the thrust axis direction was required to satisfy [ cos#t| < 0.9. 
The complete event selection has an efficiency of about 88%, and selected a total of 2 390 221 events 
in the 1993-95 data and 754 372 events in the 1991-92 data. 

Charged tracks and electromagnetic calorimeter clusters with no associated track were combined 



into jets using a cone algorithm [18], with a cone half angle of 0.65 rad and a minimum jet energy of 
5 GeV. The transverse momentum pt of each track was defined relative to the axis of the jet containing 
it, where the jet axis was calculated including the momentum of the track. 



Monte Carlo simulated events were generated using JETSET 7.4 [19] with parameters tuned by 



OPAL |2C[. The fragmentation functions of Peterson et al. [21] were used to describe the fragmentation 
of b and c quarks. The generated events were passed through a program that simulated the response 



of the OPAL detector [22] and through the same reconstruction algorithms as the data. 



2 A right handed coordinate system is used, with positive z along the e~ beam direction and x pointing towards the 
centre of the LEP ring. The polar and azimuthal angles are denoted by 9 and <f>, and the origin is taken to be the centre 
of the detector. 
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4 Tagging bb events 



Two methods were used to tag bb events, based on displaced secondary vertices and high momentum 
leptons. The first method exploits the long lifetime, hard fragmentation, high decay multiplicity and 
high mass of b hadrons. An attempt was made to reconstruct a significantly separated secondary 
vertex in each jet of the event. If a secondary vertex was found, an artificial neural network was used 
to further separate b decays from charm and light quark background. The neural network has five 
inputs, derived from decay length, vertex multiplicity and invariant mass information. The algorithm 
is described fully in p3[. For this analysis, a hemisphere was defined to be tagged if any jet in 



the hemisphere had a secondary vertex with tag variable B > 1.2 [23|, giving a hemisphere tagging 
efficiency of 32 % in bb events, and a non-b impurity of 12 %. 

Electrons and muons with momentum p > 2 GeV and transverse momentum pt > 1 GeV were 
also used to tag bb events^. Electrons were identified in the polar angle region |cos#| < 0.9 using 



a re-optimised version of the neural network algorithm described in [24|. The identification relies 
on ionisation energy loss (dE/dx) measured in the tracking chamber, spatial and energy- momentum 
(E/p) matching between tracking and calorimetry, and the multiplicity measured in the presampler 



detectors. Photon conversions were rejected using another neural network algorithm [24]. Muons were 
identified in the polar angle region | cos#| < 0.9 by requiring a spatial match between a track recon- 
structed in the tracking detectors and a track segment reconstructed in the external muon chambers, 



as m 



The tagged lepton hemispheres were further enhanced in semileptonic b decays by using informa- 
tion from the lepton p and pt, and its degree of isolation from the rest of the jet, in a neural network 



algorithm [26|. This suppresses contributions from cascade (b — > c — > £) decays (which have the wrong 
correlation of lepton sign with b hadron production flavour), charm and fake lepton background. The 
output S of the network was required to be greater than 0.7, giving a lepton sample about 75 % pure 
in b — > I decays. The lepton tags contribute an additional 5 % to the hemisphere b-tagging efficiency, 
bringing the total to about 37 % with an impurity of 13 % in Monte Carlo. 

The events used for the final analysis are those tagged by either of the b tagging methods de- 
scribed above (referred to as the T-tag) in one hemisphere, and passing the measurement selection 
requirements in the other hemisphere. These latter requirements are described below in Section rO 



for the b hadron lifetime analysis and Section 7.1 for the CP(T) analysis, and are referred to as the 



M-tag. Both hemispheres are used as measurement hemispheres in events tagged by both tags in both 
hemispheres. 

The bb purity of the combined tag T-M samples (tagged by the T-tag in one hemisphere and the 
M-tag in the other) were determined by applying an extension of the double tagging technique used 
for measuring Ry, |23]| . The number of hemispheres Ni tagged by tag i (i=T or M), and the number 
of events Nij tagged by tag i in one hemisphere and tag j in the other hemisphere, are related to the 
total number iVhad of events passing the event selection by: 



Ni = 2iV had {^ Rb +4Rc + ef s (1 - R h - R c )}, (2) 
= (2 — Sij)N hSud {Cfj e i ej R^ + e\e c ^ R c + Cf^ s e^ ds e^ ds (1 — R^ — R c )}. 



Here gives the efficiency of tag i to tag hemispheres of flavour q (q=uds, c or b). The correlations 
Cft are defined by C?- = e 1j/( e i e ]) where e^- is the efficiency to tag a qq event simultaneously with tag 
i in one hemisphere and tag j in the other. Deviations of C?- from unity account for the fact that the 
tagging in the two hemispheres is not completely independent, there being small efficiency correlations 
for both physical and instrumental reasons. The quantities i?b and R c are the fractions of hadronic 
Z° decays to bb and cc, and were taken to be R h = 0.2170 ± 0.0009 and R c = 0.173 ± 0.005 @]. 

The values of and eJp ds (which are known to be well modelled in Monte Carlo (2j|) together 
with all the correlation terms C?-, were determined from the Monte Carlo, and the values of iVhad) 

3 The notation c = 1 is employed in this paper. 
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Data Sample 


1993-95 


1993-95 


1991 


-92 




b lifetime 


CP(T) 


CP(T) 


Number of events 


2 390 221 


2 390 221 


754 


379 


Number of T-M events 


10 532 


293 416 


100 


703 


Combined tag purity (%) 


94.8 ± 1.5 


87.9 ±2.6 


01. It 


P. 

O.D 


Statistical error 


0.7 


0.4 




7 


Systematic errors: 










T-tag uds efficiency 


0.1 


1.3 




1.7 


T-tag charm efficiency 


0.3 


0.4 




0.7 


Correlations (udsc events) 


1.2 


2.2 




3.0 


Correlations (b events) 


0.5 


0.4 




0.4 


R c value 


0.1 


0.2 




0.3 


Total systematic error 


1.3 


2.6 




3.5 


Total error (%) 


1.5 


2.6 




3.6 



Table 1: The numbers of hadronic events, selected combined tag T-M events and tag purities for 
each of the data samples. The breakdown of statistical and systematic errors for each of the purity 
values is also given. 



Ni and iVy measured from the data. The five equations |2| were then solved using a \ 2 minimisation 
procedure to give the b-tagging efficiency e T of the T-tag, and all the tagging efficiencies e^, and 
e^ s of the M-tag. The b-purity Htm of the combined tag T-M sample was then calculated as: 

■q tm Ct M £t £ m -fib ^ 

C"tM e T e M ± CtM e T e M + CtM e T S£ M S 0- ~ -fib ~~ -fic) 

The results of this procedure are given in Table |l] for the three data samples: 1993-95 data with the 
M-tag used to measure the b hadron lifetime, and 1993-95 and 1991-92 data with the M-tag used in 
the CP(T) analysis. In the last case, the T-tag used vertexing in the r-<p plane only p3], and the 
purity is correspondingly lower. The purity is somewhat higher for the b lifetime analysis, as the 
requirements on the M-tag secondary vertex also provide some rejection of non-b background. In 
contrast, the CP(T) M-tag has almost equal efficiency for all flavours. 

The systematic errors resulting from each of the inputs used in the fits for the b purity are also 
given m Table |. They were evaluated using the methods and parameter ranges discussed in []23| . The 



efficiency errors include the effects of charm and light quark physics uncertainties, tracking resolution 
and lepton identification. The hemisphere tagging correlations for bb events are slightly larger than 
those in p3| , and a systematic uncertainty of ±0.02 on the C^- values is estimated, in addition to 
the Monte Carlo statistical errors. The correlations are larger because of stronger geometrical effects 
at large values of | cos 6t \ ■ The uncertainties on the uds and cc correlations C^ s and Cfj include 
systematic errors of ±0.2 but are still dominated by Monte Carlo statistics. The error from uncertainty 
in the value of i?b is negligible. 



5 Tagging the b production flavour 

The production flavour (b or b) of the b hadron in the T-tagged hemisphere was determined using 
up to three pieces of information in each event: the momentum-weighted average track charge, or 'jet 
charge'; the charge of a secondary vertex reconstructed in the hemisphere; and the charge of a high 
momentum lepton in the hemisphere. The jet charge can be calculated for every tagged hemisphere, 
and either or both of the vertex or lepton charges is always available as a consequence of the b-tagging 
techniques used. The three charges were combined using a neural network algorithm to produce a 
single production flavour tag variable Qt for each T-tagged hemisphere. 
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The jet charge Qj e t was calculated for the highest energy jet in the T-tag hemisphere as: 



where p\ is the longitudinal momentum component with respect to the jet axis and qi the charge (±1) 
of track i, and the sum was taken over all the tracks in the jet. The parameter k was set to 0.5 to 
optimise the separation between hemispheres containing b and b hadrons, including the effects of B° 



and B s mixing [27]. 



For hemispheres tagged by a secondary vertex, the charge of this vertex Q v tx was calculated as: 

Qvtx = Wi 1 i ( 5 ) 



and the uncertainty <7Q vtx as: 



Vvtx 



J2wi(l-Wi)qf (6) 



where u>, is the weight for track % to have come from the secondary, rather than the primary, ver- 
tex El. The weights Wi were obtained from a neural network algorithm using as input the track 



momentum, transverse momentum with respect to the jet axis, and impact parameters with respect 
to the reconstructed primary and secondary vertices, as in [27]. A well reconstructed vertex charge 



(with small CQ vtx ) close to +1 (—1) indicates a B + (B~), tagging the hemisphere as containing a b 
(b) quark, whilst a vertex charge close to zero indicates a neutral b hadron (e.g. B° or B°), giving no 
information on the production flavour. A vertex charge with large <JQ vtx cannot distinguish between 
Qvtx = or Qvtx = ±1) so again provides no information on the production flavour. 

For hemispheres tagged by a lepton, the lepton charge gives an estimate of the b production 
flavour, diluted by B° and B s mixing, cascade b decays, charm decays and fake leptons. The product 
QiS of the lepton charge Qi and the output S of the neural network used to select b — ► £ decays (see 
Section |j) was used as a tagging variable analogous to Qj e t and Qvtx- Higher values of S indicate 
greater probability of the lepton originating from a semileptonic b decay. 

The available production flavour estimates for each T-tagged hemisphere were combined into a 
single estimate Qt using a neural network algorithm, as in |27fl . The neural network has four inputs: 
the jet charge Qj e t> the vertex charge Q v t x and its error o"g vtx , and the lepton tag QiS. Separate 
neural networks with fewer inputs were trained for use in hemispheres with only a vertex or a lepton. 
The variable Qt is defined such that: 

_ N h (x) - N 6 (x) 
VT N h (x) + N^x) 

where N^(x) and N^(x) are the numbers of Monte Carlo b hadron and b hadron hemispheres with a 
particular value of the neural network output x. Thus hemispheres with Qt = +1 are tagged with 
complete confidence as b hadrons, hemispheres with Qt = — 1 with complete confidence as b hadrons, 
and hemispheres with Qt = are equally likely to be b or b hadrons. The modulus |Qt| satisfies 
|Qt| = 1 — 2r), where r\ is the 'mis-tag probability', i.e. the probability to tag the production flavour 
incorrectly. 

For the CP(T) violation analysis, both the production flavour estimate from the T-tag hemisphere 
and the jet charge Qj e t in the M-tag hemisphere are used to infer the production flavour of the b 
hadron in the M-tag hemisphere. In this hemisphere, the parameter k in equation |I| is set to zero, so 
the jet charge becomes simply the average of the charges of the tracks in the jet. This avoids being 
sensitive to the decay flavour of mixed or unmixed B° and B s mesons, but is still sensitive to the 
production flavour of the b hadron via the information carried by the fragmentation tracks in the 
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jet p9| . This jet charge is used to produce a tagging variable Qm defined similarly to Qt for the T-tag 
hemisphere. The information from both hemispheres is combined into a single variable Q2, defined as: 



Q-2 



(1-Q T )(1 + Q 



U) 



(1-Qt)(1 + 0m) + (1 + Qt)(1-Q 



The variable Q2 is positive (negative) if the event is tagged as containing a b (b) hadron in the M-tag 
hemisphere. 

The Qjet and Q v tx distributions are not charge symmetric because of detector effects causing a 
difference in the rate and reconstruction of positive and negative tracks. These effects are caused by 
the detector material and the Lorentz angle in the tracking chambers |pLCf| . They were removed by 
subtracting offsets from the Qj e t and Q v tx values before the calculation of Qt and Qm- The offsets 
were determined from data using the inclusive T-tagged samples, and were found to be fractions 
0.018 ± 0.002, 0.034 ± 0.002 and 0.028 ± 0.002 of the RMS width of the Q jet (rt = 0.5), Q jet (« = 0) and 
Qvtx distributions, respectively. 

The distributions of Qt and Qm in T-tagged data and Monte Carlo events are shown in Figures |l|(a) 
and 0(b). The sharp drop at about \Qt\ = 0.8 in Figure [l](a) is due to the irreducible fraction of 
lepton tagged events that are tagged incorrectly due to B° and B s mixing. Around 30 % of events 
have equal numbers of positive and negative tracks in the jet in the M-tag hemisphere, giving zero jet 
charge and Qm = 0. These events are not shown in Figure ||(b). Some discrepancies between data 
and Monte Carlo distributions are visible. These discrepancies are not important, provided that the 
Qt values in data correctly describe the data mistag probabilities. To check that this is the case, 
events where both hemispheres are T-tagged, yielding b production flavour estimates Qti and Qt2, 
were used. Since the two hemispheres must contain b hadrons produced with opposite flavour, the 
product £tt = QtiQt2 is negative if both (or neither) hemispheres are tagged correctly, and positive 
if only one is tagged incorrectly The distribution /(£tt) of £tt in data thus allows the production 
flavour estimate Qt to be checked. The function: 

represents the 'wrong sign fraction' at a particular value of |£tt|, and should satisfy <7(£tt) = (1 — 
£tt)/2 if Qt correctly describes the mis-tag probabilities. The distribution 5(£tt) in data is shown in 
Figure 0(c), together with a linear fit. The distribution has the expected form, and the fitted gradient 
is dg/d^TT = —0.514 ± 0.012, showing that the average magnitude of Qt is correct to a relative 
precision of 2.6 %. 

The same technique was used to study Qm, via the 'cross-tag' product £tm = QtiQm2, as Qm has 
a smaller tagging power than Qt- The function <?(£tm) was defined in an analogous way to <?(£tt) in 
equation 0. The distribution <?(£tm) in data is shown in Figure 0(d), and the corresponding gradient is 
—0.545 ± 0.018. Since Qt has already been shown above to be correct to 2.6 %, all of this discrepancy 
is attributed to Qm, which therefore has a bias of up to about 10 %. These values are used to calculate 
tagging systematic errors. 



6 B + and B lifetime analysis 

The bb event tagging and production flavour tagging described above is common to both the b hadron 
lifetime and CP(T) analyses. The remainder of each analysis — the M-tag and the fits to the data — are 
specific to each analysis. The b hadron lifetime analysis is described in this section, and the CP(T) 
analysis in Section 0. 
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Figure 1: Production flavour tagging of bb events: (a) Distribution of the flavour tagging variable 
Qt in data (points) and Monte Carlo (histogram) T-tagged hemispheres. The contributions from b 
hadrons, b hadrons and non-b backgrounds are shown, (b) The analogous distributions for the M-tag 
hemispheres. The 30 % of the sample (90000 events in data) with Qm = are not shown. In both (a) 
and (b), the error bars on the data points are smaller than the symbols, (c) Distribution of the wrong 
sign fraction <?(£tt) in data double T-tagged events, with a linear fit superimposed, (d) Distribution 
of the wrong sign fraction ^(^tm) in data double T-tagged events, again with a linear fit superimposed. 
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6.1 b hadron reconstruction 



The b hadron reconstruction used for the M-tag in this analysis is similar to that used in ||. It 
aims to reconstruct a relatively small sample of clear secondary vertices where each track can be 
unambiguously associated to either the primary or secondary vertex. 

The algorithm considers all tracks in a jet which have momentum p > 0.5 GeV, impact parameter 
(in the r-<p plane) | cZo | < 1cm and error on the impact parameter ad Q < 0.1cm. All possible sets of 
assignments of these tracks to the primary and secondary vertices ('arrangements') were considered, 
requiring at least two tracks to be assigned to the secondary vertex, and including the combination 
where no tracks at all are assigned to the primary vertex. 

The positions of the two vertices in each arrangement were determined by fitting all the assigned 
tracks to a common vertex in three dimensions. For the primary vertex, an additional constraint from 
the average beam spot position was used, determined from a fit to the tracks in many consecutive 
events p0j |. The x 2 °f the arrangement was calculated as the sum of the x 2 values for the primary 
and secondary vertex fits, and the fit probability of the arrangement determined from the x 2 value 
and the number of degrees of freedom in the two vertex fits. 

For the jet to be accepted as having a clear secondary vertex, the following conditions had to be 
satisfied: 

1. One and only one arrangement (the 'best arrangement') has a fit probability exceeding 1 %. 

2. All other arrangements have a x 2 value exceeding that of the best arrangement by at least 4. 

3. The decay length L of the secondary vertex, divided by its error cr^, satisfies L/ai > 3. The 
decay length is calculated from the distance between the fitted primary and secondary vertices, 



using the direction of the jet axis as a constraint [23| 



4. The decay length must be positive, but less than 3 cm. The decay length is positive if the 
secondary vertex is displaced from the primary vertex in the same direction as the jet momentum 
vector. 

These requirements ensure that the best arrangement has an acceptable probability, that no other 
arrangement is likely to be the correct one, and that the primary and secondary vertices are well 
separated. Jets containing b hadrons with short decay lengths will tend to fail requirements 2 or 3, 
whilst those with a mis-measured track which does not fit with either the primary or secondary vertex 
will tend to fail requirement 1. 

The charge Q of the secondary vertex and its error gq were then calculated, using equations ||| 
and ^, with track weights Wi optimised for this vertex finding algorithm. To ensure the vertex charge 
was well reconstructed, the error <jq was required to be less than 0.7. All tracks in the jet were used 
in calculating the vertex charge, including those which failed the tighter selection used for the initial 
vertex finding. 

A total of 10 532 combined T-M tagged events were found in the data with this selection for the 
M-tag. The distribution of the M-tag vertex charge, together with the Monte Carlo prediction, is 
shown in Figure B Clear peaks at Q = ±1 and are seen, corresponding to high purities of charged 
and neutral b hadrons. In the Monte Carlo, the neutral b hadrons are a mixture of approximately 
63% B°, 24% B s and 13% b baryons (denoted by Ab). The charged b sample is almost entirely B + , 
with a very small contribution from charged b baryons (H^) which is estimated at about 1 % and is 
neglected. 

6.2 Excess proper time reconstruction 

The decay length distributions of secondary vertices selected by the above algorithm in data and 
Monte Carlo are shown in Figure ||(a). The decay lengths are very biased towards large values, since 
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Figure 2: Distributions of the vertex charge Q in the data (points) and Monte Carlo (histogram), 
normalised to the same number of entries. The expected contributions from charged b hadrons, 
neutral b hadrons and non-b background are indicated. 
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the algorithm preferentially selects well separated secondary vertices. Therefore in order to extract 
the b hadron lifetimes, the excess decay length method is used, as in |E, 13 1. For each selected event, 



the minimum b hadron decay length that would still result in a resolvable secondary vertex passing 
all the requirements described in Section ^O] was determined. To find this minimum decay length, 
all the tracks assigned to the secondary vertex in the best arrangement were assumed to come from 
the b hadron. They were then translated along the direction of the jet axis towards the primary 
vertex (as if the b hadron had decayed at an earlier time), and the entire vertex selection repeated 
(requirements 1-4 of Section 6.1 and oq < 0.7, including recalculating the x 2 values for all possible 
vertex arrangements). The translation distance at the point where the modified arrangement just fails 
one of the requirements defines the excess decay length L ex . The distributions of L cx in data and 
Monte Carlo are shown in Figure §(b). 

For most selected vertices, the point of failure occurs when one of the less good arrangements has 
an improved x 2 , and fails requirements 1 or 2. Others fail because the L/a separation becomes too 
small or the vertex charge error oq becomes too large. 

The distribution of excess decay length is approximately a negative exponential, with the same 
slope as the b hadron decay length distribution before the vertex selection requirements were imposed. 
The distribution is exactly an exponential provided that the assignment of b hadron decay tracks to 
the secondary vertex and fragmentation tracks to the primary vertex is correct, and that the effects 
of vertex resolution and the lifetime of the charm hadron produced in the b hadron decay can be 
neglected. 

That the distribution is exponential is most easily seen in terms of the excess decay proper time, 
obtained from the excess decay length via the reconstruction of the b hadron energy. The effect of 
the charm hadron lifetime will be considered first. Ignoring resolution effects, the rate of events F(t) 
with excess decay time t is given by the convolution of the lifetime exponentials e~ tb ^ Tb and e~ tc ^ Tc 
for the decaying b and charm hadrons. Here i& is the decay time and the lifetime of the b hadron, 
and similarly for the charm hadron. The convolution is obtained by integrating over the excess decay 
time t' of the b hadron, defined as t' = tf, — to, where to is the minimum decay time for this event to 
pass the vertex selection requirements. Since the introduction of excess decay time just corresponds 
to the redefinition of zero time at t& = to and does not affect the form of the b lifetime exponential, 
the distribution F{t) is given by replacing % by t', t c with t — t', and integrating: 



F(t) oc /' 

J a 



* e -t'/T be -(t-t')/Tc dt > 



oc e- t/Tc 



e - t '(i/ T{ _l/ Tc )1*- 6 

J t'=a 



where normalisation factors have been neglected. The upper limit b is simply given by b = t, as 
only b hadrons decaying with excess time smaller than t can contribute. The lower limit a = —to, 
corresponding to b hadrons decaying at tf, = 0. As long as this minimum b hadron excess decay time 
is in magnitude much larger than the maximum contributing charm decay time, the lower limit can be 
approximated by a = — oo. The requirements of Section [O] select events with a single secondary vertex 
well separated from the primary vertex. This ensures a long b hadron decay time, and suppresses 
events with a resolvable tertiary vertex from a long-lived charm hadron, effectively truncating the 
charm decay exponential. These conditions ensure that the approximation a = —oo is valid, and the 
integral finally becomes 

F(t) oc e~ t/Tb . 

A similar argument holds for the effect of the finite detector resolution. In this case, the resolution 
function replaces the charm hadron lifetime exponential in the convolution, and the limits become 
a = — oo and b = oo. Since the convolution of an exponential with any finite function is another 
exponential with the same decay constant as the original, the b hadron lifetime distribution is again 
recovered. 
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Figure 3: Distributions of (a) decay length L and (b) excess decay length L cx in data (points with error 
bars) and Monte Carlo (histogram), for all selected secondary vertices. The expected contributions of 
b hadrons and non-b background are indicated. 
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However, for a small fraction of events (mainly those with low L eK ), one or more tracks are mis- 
assigned from the primary to the secondary vertex or vice versa, which introduces distortions in the 
excess decay length and proper time distributions. Such events also tend to have some tracks with 
vertex charge weights u;, ~ 0.5 (i.e. not clearly assigned to either vertex), and so are concentrated in 
the regions of Q away from the peaks at integer values. 

The excess decay length was combined with an estimate of the b hadron energy E^ in each event 
to determine the excess proper time t, via the relation: 

_ m h L cx 

where m\, is the b hadron mass, taken to be that of the B+ and B°, i.e. 5.279 GeV @. The b hadron 
energy was estimated using a technique described in [^T|. First, the energy of the b jet £?bjet was 
estimated by treating the event as a two-body decay of a Z° into a b jet of mass and another 
object making up the rest of the event. Then, the charged and neutral fragmentation energy -Ebfrag 
in the b jet was estimated using the charged track weights Wi, and the unassociated electromagnetic 
calorimeter clusters weighted according to their angle with respect to the jet axis. Finally the b hadron 
energy was calculated as = -Ebjet — ^bfrag- 

The resulting distributions of b hadron energy in data and Monte Carlo are shown in Figure ^a). 
The agreement is generally good, and the small differences around the peak are within the uncertainties 
due to the imprecise knowledge of b fragmentation. In addition to all the vertex requirements described 
above, selected M-tag hemispheres were required to have E^ > 20 GeV, which suppresses to a negligible 
level events where the M-tagged jet originates from a gluon and not a b quark. This cut is included 
in the definition of the M-tag and is included in the event counts and purities of Table ||. The energy 
resolution in Monte Carlo bb events is shown in Figure ||(b). The reconstructed energy has a mean 
equal to the true energy, but has asymmetric tails. However, these tails are small enough not to have 
a significant effect on this analysis. 

6.3 Fit and results 

The lifetimes of the B+ and B° mesons were extracted by using a maximum likelihood fit to the mean 
excess proper time t as a function of the modulus of the vertex charge Q. All events with < t < 15 ps 
were used in the fit. The data were divided into ten bins between \Q\ = and \Q\ = 2, and all events 
with \Q\ > 2 put into an eleventh bin. The mean excess decay time (t)j was then calculated for each 
\Q\ bin j, and compared to the fit prediction tj. The latter depends on the lifetime t s of each source 
s (b hadron type or background) and the fraction of each source /| expected in bin j. 

The fractions /| depend on the vertex charge Q. If the charge tagging were perfect, only charged 
b hadrons would be expected at Q = ±1, and only neutral b hadrons at Q = 0. However, there is 
some cross-contamination, as can be seen in Figure As the measurement of the B + and B° lifetimes 
depends crucially on the level of this contamination, this was fitted from the data by exploiting the 
correlation of the vertex charge Q with the opposite (T-tag) hemisphere b production flavour estimate 
Qt- For example, a b hadron in the T-tag hemisphere implies a b hadron^ (B~, B°, B s or Ab) in the 
M-tag hemisphere, giving a correctly reconstructed charge of Q = — 1 or Q = 0. The number of such 
events reconstructed with Q = +1 therefore gives information on the number of neutral b hadrons 
incorrectly reconstructed as charged b hadrons, since a true B + (being a b rather than a b hadron) 
cannot be opposite another b hadron. In fact, the charged/neutral separation is a function of the 
continuous variable Q, and is described by a parameterisation constrained by the above correlation. 
These parameters were determined at the same time as the source lifetimes t s as additional parameters 
in the fit. 

4 By convention, the b mesons considered as particles (B + , B°, B s ) contain a b antiquark, and the antiparticles (B~, 
B°, B s ) contain a b quark. The opposite is true for baryons, so a Ab contains a b quark and a Ab contains a b antiquark. 
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Figure 4: (a) Reconstructed b hadron energy distributions for data (points with error bars) and Monte 
Carlo (histogram). The position of the minimum energy cut is shown by the arrow, (b) Resolution of 
the reconstructed b hadron energy in Monte Carlo bb events, after the minimum energy cut has been 
applied. 
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In more detail, the total likelihood of the event sample was given by: 

£ = £ time.-Q £ tag (§) 

% 

where £ time represents the likelihood from the fit to the lifetime as a function of vertex charge Q and 
£* ag represents the likelihood of each event used to determine the charged/neutral separation. The 
former is calculated in bins of Q, whilst the latter is determined event by event, and the product is 
taken over all events i. 

The logarithm of the time likelihood £ time is given by: 

m£ ti m e = E _ iV ^(|i +lnr ^ 

where the index j runs over the bins of vertex charge Q. The term inside the sum represents the 
log-likelihood to measure a mean decay time (t)j in a sample of Nj events distributed according to 
a negative exponential with lifetime Tj. Although the data events are not distributed according to a 
single exponential, the differences are sufficiently small that this expression can be used. The expected 
true mean decay time in bin j is given by: 

T j = D((Q) j ) + Y,f J sT s (9) 

s 

where fl is the fraction of events from source s in bin j, and t s is the lifetime of source s. There are 
nine sources in total: 

s = {B+,B-,B°,B°,B s ,B s ,A b ,A b , background} 

and the lifetimes of particle and antiparticle are assumed to be equal (r B + = r B - etc). The background 
is characterised by a lifetime Tb g , taken from Monte Carlo. 

The function D((Q)j) in equation || accounts for the distortions in the excess decay length distri- 
bution caused by mis-assigned tracks discussed in Section |6.2| . These can be seen in Figure [|, which 
shows the mean excess proper time (t) as a function of \Q\ for both data (points) and Monte Carlo 
(solid line). Events which have one or more mis-assigned tracks tend to have smaller than average 
excess proper time, and to be concentrated away from integer values of Q. This reduces the mean 
excess decay time in these regions, and consequently increases the mean time close to integer values of 
Q. This effect, seen clearly in both data and Monte Carlo, is parameterised by the periodic correction 
function D, which has the form: 

D(q) = d(q - 1/2) 2 -d/8 {for < q < 1} (10) 

where q = {Q — int(Q)) is the non-integer part of Q. The parameter d characterises the amplitude of 
the distortion, and is left as a free parameter in the fit. The functional form of this correction was 
chosen by studying the effect in Monte Carlo. 

The second part of the overall likelihood in equation || is the tag likelihood £ tag . It is given for 
each event i by: 

s 

The function P®{Q) gives the probability of each source s as a function of the b hadron vertex charge 
Q. For the B + this function is given by: 

f (l-/bg)c 3 forQ<-l 

P Q (n\ = ) (l-/bg)(co + (c -c 3 )Q) for-KQ<0 

B+W) I (l-/bg)(co + (ci-c -c 2 )Q + C2 Q 2 ) forO<Q<l 

{ (1 - A g )ci for Q > 1 
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Figure 5: The distribution of mean excess decay time (t) as a function of the modulus of the vertex 
charge \Q\. The data are shown by the points with error bars, and the Monte Carlo (reweighted to 
the same B + and B° lifetimes as measured in the data) distribution is shown by the solid line. The 
prediction of the fit is shown by the dashed line. The bin between \Q\ = 2 and \Q\ = 2.2 contains all 
events with \Q\ > 2. 



where /b g is the fraction of non-bb background, cq represents the fraction of signal B + events at 
Q = 0, and c\ the fraction of signal B + events at Q = 1. A quadratic interpolation with coefficient C2 
is used in the range < Q < 1, and the fraction of B + is constant for Q > 1. A small fraction C3 of 
B + events has Q < —1, and a linear interpolation is used for the fraction of B + between —1 and 0. 
This functional form is illustrated in Figure ^ and is seen to give a reasonable description of the B + 
fraction in Monte Carlo. Large variations in the parameters cq and C3 are considered when assessing 
the systematic errors. 

The fraction of B~ is given by charge symmetry: P®_(Q) = P® + (—Q). The remaining signal 
fraction is neutral b hadrons, so the fractions of B° and B° are given by: 

p b° (Q) = p£o(Q) = (i - fl - flja - /b g - p§AQ) - p§-(Q)) (12) 

where /g s and are the fractions of B s and At, in the neutral b hadron sample. The values for 
B s and Ab are given in an analogous way to equation |l2| but with the factor (1 — /g s — /? ) replaced by 

/b s or /A b - Finally the background function P® g (Q) = /b g , i.e. the background fraction is a constant 
independent of Q, as found in Monte Carlo. 

These functions together describe the source fractions as a function of Q in terms of four param- 
eters: ci and C2, which are left free in the fit; and co and C3, which are input from Monte Carlo. 
The values of c\ and C2 are constrained by the correlation between the type of b hadron in the M-tag 
hemisphere and the production flavour (b or b) of the other b hadron in the T-tag hemisphere. The 
contamination of neutral b hadrons at Q ~ ±1 (which is given by 1 — c\ — C3) is thus determined 
almost entirely from the data, whilst the contamination of charged b hadrons at Q ~ (given by Co) 
has to be taken from Monte Carlo, since the vertex charge provides no distinguishing power between 
neutral b and b hadrons. 
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Figure 6: Fraction of B + vertices as a function of vertex charge Q in Monte Carlo bb events (points 
with error bars). The parameterisation is shown by the solid line, and the levels Co, c\ and C3 by the 
dotted lines. The quadratic coefficient C2 is used between Q = and Q = 1. Note that the non-bb 
background fraction is not included, and the B + fraction is plotted using a logarithmic scale. 



The function PJ^Qt) in equation [Tl] gives the probability for each source to be tagged by the 
opposite hemisphere flavour tag of value Qt: 



(1-Q t )/2 for s = B+,B°,B s ,A t 
(1 + Qt)/2 for s = B-,B°,B s ,A t 
1 for s =background 



Finally, the fractions /| of each source contributing in each bin j of the time likelihood (equation ||) 
are given by the average of the source fractions /j for each of the events i in bin j, where: 

4 _ p?(Qi)*T(.QTi) 

Js 



Y, s >P3{Qi)Pj>{QTi) 

In total, 5 parameters were left free in the fit: the B + and B° lifetimes t-q+ and Tgo, the charge 
separation parameters c\ and C2, and the distortion parameter d. The B s and Ab lifetimes were taken 
to be 1.54 ± 0.07 ps and 1.22 ± 0.05 ps |Q. The background lifetime was taken to be that in the Monte 
Carlo, Tb g = 2.1 ps, and the background fraction was measured to be /b g = 0.052 ±0.015 (see Table |l]). 
The remaining parameters /g s , /° , cq and C3 were taken from the Monte Carlo. Uncertainties on all 
these parameters are considered when evaluating the systematic errors. 

The results of the fit for the B + and B° lifetimes are: 

r B + = 1.643 ± 0.037 ps 
r B o = 1.523 ± 0.057 ps 

where the errors are statistical only. The correlation coefficient between the two measured lifetimes is 
—0.53, and their ratio is TB+/ r B° = 1-079 ± 0.064. The other fitted parameters were measured to be 
ci = 0.84 ± 0.03, c 2 = 0.10 ± 0.18 and d = 1.76 ± 0.21 ps. 
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Source 


Atr+ (ps) 


Atro (ps) 


Atr+ / TrO 


B s lifetime 


0.000 


0.026 


0.019 


b baryon lifetime 


0.000 


0.011 


0.008 


B s /B° fraction 


0.000 


0.001 


0.001 


Ab/B° fraction 


0.000 


0.025 


0.018 


Background lifetime 


0.005 


0.009 


0.003 


Background fraction 


0.010 


0.015 


0.004 


Fit procedure 


0.014 


0.023 


0.022 


Charge separation 


0.005 


0.019 


0.015 


b fragmentation 


0.008 


0.009 


0.000 


Detector resolution 


0.008 


0.007 


0.009 


Silicon alignment 


0.009 


0.011 


0.012 


Total 


0.025 


0.053 


0.041 



Table 2: Systematic errors on the measured values of r B +, r B o and t b +/t b o. The values for t b +/t b o 
take into account correlations between r B + and r B o . 



The distribution of mean excess decay time (t) as a function of |Q| is shown in Figure [|, for 
data and for Monte Carlo simulation reweighted to the measured B + and B° lifetimes. The Monte 
Carlo sample is 13 times larger than the data sample, and the x 2 between data and Monte Carlo 
distributions is 9.4 for 8 degrees of freedom. The results of the fit to the data (the values of Tj) are 
shown by the dashed line, and the fit follows the data reasonably well. The distribution of data excess 
proper time in various regions of \Q\ is shown in Figure ^, again together with the reweighted Monte 
Carlo simulation. The deviations from a pure exponential away from integer values of Q can clearly 
be seen in both data and Monte Carlo. 



6.4 Systematic errors 

The systematic errors on the B + and B° lifetime measurements and their ratio are summarised in 
Table |2], and discussed in more detail below. In most cases, the systematic errors are evaluated by 
repeating the fit in either data or Monte Carlo, varying the appropriate parameter. 

B s and Ab contamination: The neutral b hadron sample consists not only of B° mesons but an 
irreducible contribution from B s and Ab- Uncertainties in both the lifetimes (r Bs and T\ b ) and 
sizes (/ Bs and /° ) of these backgrounds affect the B° lifetime measurement, the B s lifetime 
being the largest single source of uncertainty. The efficiencies for B°, B s and A^ to pass the 
vertex selection requirements are similar, but somewhat sensitive to their average charged decay 
multiplicities. However, by far the largest uncertainty in the level of contamination comes 
from the knowledge of the B s and Ab production fractions in Z° — > bb events, taken to be 
/(b B s ) = 10.5j£? % and /(b A b ) = 10.1+1;? % Q. 

Non-bb background contamination: The effective mean decay time of the non-bb contamination 
in the selected sample is 2.1 ps in Monte Carlo. This has two distinct components — a contribution 
from mis-measured tracks and strange particle decays causing fake b decay vertices with long 
decay lengths, and a contribution from the decay of charm hadrons which have genuine lifetime. 
Charm events have both contributions, giving an effective lifetime of about 1.8 ps. Only the 
former contribution is present in uds events, which have an effective lifetime of about 2.3 ps. 

To study the modelling of the background lifetime, two control samples were used, generated 
by modifying the selection for the T-tag. Requiring the vertex tag variable B in the T-tag to 
be between and 1.2 gave a sample consisting of 9% uds, 18% cc and 72% bb events, and 
requiring the T-tag hemisphere to fail the vertex tag pre-selection [23 gave a sample of 41% 
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Figure 7: Distributions of excess proper time in four regions of absolute vertex charge \Q\, for data 
(points with error bars) and Monte Carlo simulation (histogram). The Monte Carlo B+ and B° 
lifetimes have been reweighted to the same values as measured in the data. 
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uds, 23 % cc and 35 % bb events. The effective lifetimes measured in the M-tag hemispheres for 
these samples were found to agree in data and Monte Carlo to better than 0.05 ps. From this, 
an upper limit due to the mis-modelling of the background in the primary b-tagged sample was 
estimated as 0.1 ps, and used to assess the systematic error on the B + and B° lifetimes. 

Fit procedure: The entire fitting procedure, including the derivation of the excess decay lengths and 



the correction for the distortion given in equation 10, was tested on a fully simulated Monte Carlo 
sample 13 times larger than the data sample. In this sample, where the B + and B° lifetimes 
were both 1.6 ps, the fit gave the results t b + = 1.588 ± 0.014 ps and r B o = 1.623 ± 0.017 ps. The 
larger of the deviations of these results from the true values or the statistical errors were taken 
as systematic errors due to the fitting procedure. As can be seen from Figures [| and |?], the 
Monte Carlo provides a good description of the time distributions in the data. 

Additional Monte Carlo studies were performed to check the correctness of the fit procedure 
and errors. To verify the errors returned by the fit, it was performed on many Monte Carlo 
subsamples, and the distribution of fitted lifetimes studied. The Monte Carlo was also reweighted 
to change the B + and B° lifetimes one at a time. In each case, the fit correctly recovered the 
modified lifetime, whilst returning a stable result for the lifetime of the unmodified b hadron. 

As a cross check, the fit to the data was repeated with the parameterisation of the distortion 
given in equation |l^ replaced with the distortions measured in the large Monte Carlo sample. 
Consistent results were obtained, with the values of r B + and r B o changing by 0.017 ps and 
—0.013 ps respectively. 

Charge separation: The fraction of neutral b hadron contamination at |Q| ~ 1 (1 e. \ — c\ — C3) was 
fitted from the data. However, the corresponding contamination of charged b hadrons at Q ~ 
{i.e. Co), and the fraction of wrong sign charged b hadrons C3 were taken from the Monte Carlo, 
and varied by ±50% to assess the systematic error (see Figure ||). This variation is larger than 
the difference seen between data and Monte Carlo in the neutral contamination of the charged b 
sample, and also larger than the variation resulting from changing the b hadron charged decay 
multiplicities within their experimental uncertainties p2| . 

The fit to determine the charge separation depends on the T-tag hemisphere production flavour 
tag (5t> which is subject to mis-tag and jet and vertex charge offset subtraction uncertainties, as 
discussed in Section || The effect of these uncertainties on the b hadron lifetime measurement 
is negligible. 

b fragmentation: The effect of uncertainties in the average b hadron energy (x^) — -E\)/-^beam was 
assessed in Monte Carlo by reweighting so as to vary (xe) in the range 0.702 ± 0.008 [j32]| , 
and repeating the lifetime fit. The effect on the lifetimes is small since the b hadron energy is 
estimated event by event. 

Detector resolution and alignment: The error due to uncertainty in the tracking detector reso- 
lution was assessed in Monte Carlo by applying a global 10 % degradation to the resolution of 
all tracks, independently in the r-cp and r-z planes, as in p3[] . The lifetime measurements are 
also sensitive to the effective radial positions of the silicon detectors (both the positions of the 
detectors themselves and the positions of the charge collection regions within them). These are 
known to a precision of ±20 /xm from studies of cosmic ray events [Q. The resulting uncertainty 
was calculated by applying 20 /xm radial shifts to one or both silicon layers in Monte Carlo and 
repeating the lifetime fits. 

The total systematic errors amount to ±0.025 ps on the B + and ±0.053 ps on the B° lifetimes, and 
±0.041 on their ratio, where correlated systematic errors have been taken into account. 
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7 CP(T) violation analysis 



In the Standard Model framework for describing CP violation, according to the formalism given in 
equation |], an asymmetry is predicted in the time dependent inclusive decay rates of B° and B° 
mesons: 

, _ r(B° — > anything) — r(B° — > anything) 
^ ' ~ r(B° anything) + r(B° anything) 

For a totally unbiased selection of B° and B° mesons, the form of this asymmetry as a function of true 
proper decay time t' is given by: 

A(t>) = a cp {^pisin(Am d tO - ^ (^)} (13) 

where a cp is the CP-violating observable, Amj is the B° oscillation frequency and r^o the B° lifetime 
jl2,^3|. The parameter a cp is related to the CP-violation parameter ee by: 

Re(e B ) = °f (14) 

using the convention that |ee| is small. Searching for an asymmetry of the form given in equation pl| 
therefore provides a method of probing the value of e-Q. In the Standard Model, Re(eB) is expected to 
be around 10 -3 |12]| , but it could be up to an order of magnitude larger in superweak models pf| . 

In this analysis, the T-tag is used to tag bb events, and an inclusive reconstruction algorithm 
is used to reconstruct the b hadron decay time t in the M-tag hemisphere. The production flavour 
(b or b) of this M-tagged b hadron is given by the combined tagging variable Q2 (see section ||), 
using information from both hemispheres of the event. No attempt is made to separate B° and B° 
decays from other b hadron decays in the M-tag hemisphere. A similar asymmetry to that given in 



equation |13j is expected for B s mesons, with Amj replaced with the B s oscillation frequency Am s , 
but it is expected to be at least an order of magnitude smaller and is neglected. No asymmetry is 
expected for B + and At,. These other b hadrons therefore dilute the expected B° asymmetry, but do 
not change its form. Hence no attempt is made to remove B + decays by requiring well reconstructed 
neutral vertices, as this would not greatly increase the sensitivity to CP violation in the B° decays, 
and would lead to a large loss in reconstruction efficiency. 

Although the time dependent rates of B° and B° decay are predicted to be different, this does not 
violate the CPT invariance prediction of equal total decay rates. The lifetimes of B° and B° are thus 
still expected to be equal. If CPT violation were present, the lifetimes of b and b hadrons could be 
different: 

1 /At 



Tt " i 1 2 V r 

where r av is the average and (Ar/r)b the fractional difference in lifetimes. The value of (Ar/r)b 
is determined by measuring the lifetimes of b and b hadrons, using the same time reconstruction 
algorithm as for the CP-violation analysis. 

7.1 b hadron reconstruction 

The b hadron decay length was reconstructed using a secondary vertex finding algorithm similar to 



those employed in |^8| and |31]. In the highest energy jet in the M-tag hemisphere, the two tracks with 
the largest impact parameter with respect to the primary vertex were used to form a 'seed' vertex. 
All tracks consistent with this seed vertex were then added to it using an iterative procedure. Only 
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tracks satisfying the quality requirements detailed in Section 6.1 were considered for inclusion in the 
secondary vertex. 

The resulting vertex was required to have at least three tracks, and the invariant mass (assuming 
all tracks to be pions) was required to be at least 0.8 GeV. The transverse miss distance (the projection 
of the vector between primary and secondary vertices onto a plane orthogonal to the total momentum 
vector of the tracks assigned to the secondary vertex) divided by its error was required to be less 
than 3 These requirements help to eliminate badly reconstructed and fake secondary vertices. 

The decay length L between primary and secondary vertices was then calculated, using the jet axis 



direction as a constraint, as in |23j. An acceptable secondary vertex was found in approximately 70 % 
of hemispheres, for both bb and background (cc and light quark) events. 

The b hadron energy was estimated as described in section |6.2| , calculating the charged fragmenta- 
tion energy using weights Wi tuned with this alternative vertex finder. The estimate t of the b hadron 
decay proper time was then calculated as before: 

m\)L 



El -ml 



In Monte Carlo, the resolution of this estimate can be described by the sum of two Gaussians. For the 
1993-95 data, the narrower Gaussian has an RMS width of 0.33 ps and contains 65% of the events, 
and the wider Gaussian has a width of 1.3 ps. For the 1991-92 data, where the vertex reconstruction 
was done in the r-(p plane only, the narrow Gaussian has a width of 0.33 ps and contains 64 % of events, 
and the wider Gaussian has a width of 1.4 ps. Since the resolutions are similar, the two data samples 
were combined and a single resolution function used for the entire sample. These resolutions are an 
average over all true decay proper times t', and significant non-Gaussian effects can be seen in small 
slices of t' , as in [^TJ. These effects are caused by the presence of tracks from the primary vertex, and 
make a full description of the resolution as a function of t' very complicated. However, they are not 
important for the analysis described here, which does not rely heavily on an accurate description of 
the decay time resolution. 

7.2 Fit and results for €b 

The CP-violating parameter a cp was extracted using a x 2 fit to the observed time dependent asym- 
metry A{i) in 34 bins of reconstructed time in the range —2 to 15 ps. Within each time bin i, the 
observed asymmetry was calculated in ten bins j of \Q 2 \ (0 < IQ2I < 1) to make best use of the tag- 
ging information in each event. These ten estimates of the asymmetry were averaged with appropriate 
weights to calculate the overall observed asymmetry ^4° bs in each time bin i. The observed asymmetry 
was compared with the expected asymmetry A^ Tne calculated for a given a cp , taking into account the 
time resolution and dilution from non-B° decays in the event sample. 

The corrected observed asymmetry in bin i of reconstructed time t and bin j of | Q2 | is given by: 

_/V b . - /V?>. 
^obs _ iV »j iV u 



t.l 



(\Q2\)ij(Nh + N$ 



and the error cr^obs is given by: 



(7 Aohs 



({\Q2Dij A^ S ) 2 N%+N% 
2<|Q a [)« \ N*M 



where (N^) is the number of events with Q2 > (Q2 < 0). The factor l/{\Q2\)ij corrects for the 
tagging dilution (mis-tagging), which reduces the observed asymmetry for imperfectly tagged events. 
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Figure 8: Asymmetry of tagged b and b hadrons as a function of reconstructed time t in data (points) 
and fit (solid line). The expected asymmetry for a cp = 0.15 is shown by the dotted line. 



The ten estimates Af? s were then averaged, weighting according to (a A ob S ) 2 to derive the corrected 

observed asymmetry ^4° bs . 

The true asymmetry as a function of reconstructed time t is given by: 



A *»*(t) = (1 - / bg )/ B o 



/ °° P(t)A(t')R(t - t')dt' ' 
f™P(t)R(t-t')dt> 



where /b g is the fraction of non-bb events in the sample, _f B o is the fraction of B° in the bb part of the 
sample, P(t') is the lifetime exponential P(t') = (l/r B o)e _ * / r B° ; A(t') is the asymmetry as a function 
of true proper time t' given in equation |l3], and R(t — t') is the time resolution function. R(t') is the 
sum of two Gaussians, with parameters given in section |7.1| . The expected asymmetry in time bin i 
was calculated from the mean reconstructed decay time of all events in bin i: A\ rue = A true ((t)i). 

A binned x 2 fit was performed to the distribution of asymmetry as a function of reconstructed 
time. Both 1991-92 and 1993-95 data samples were used together, setting the background parameter 
to the average impurity of the two samples, /bg = 13.7 ± 3.0% (see Table Jl]). The values of r B o and 
Amd were taken to be r B o = 1.56 ± 0.04 ps and Ama = 0.464 ± 0.018 ps -1 [Q]. The parameter / B o was 
taken from the Monte Carlo to be / B o = 0.41. 

The observed asymmetry for the 394 119 data events is shown in Figure ||, together with the result 
of the one parameter fit to a cp . The fit result is 

a cp = 0.005 ± 0.055 

where the error is statistical only. 

The form of the asymmetry given in equation 13 assumes that the reconstruction efficiencies for 
all decay modes are equal. In particular, if the efficiency to reconstruct B° decays to final states 
containing no charm hadron, one charm hadron and two charm hadrons are different, then additional 
asymmetries may be seen f33]|. This is because the semi-inclusive B° decays to final states containing 
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Source 


Aa cp 


Ac cp 


R u lifptimp 

J_> ±l±t; LlIJ.lt; 




0.002 


0.000 


Arrid value 




0.001 


0.001 


B° fraction 




0.002 


0.002 


Flavour tagging 


offsets 


0.003 


0.013 


Flavour tagging 


mis-tag 


0.009 


0.005 


b fragmentation 




0.008 


0.006 


Time resolution 




0.002 


0.000 


Total 


0.013 


0.015 



Table 3: Systematic errors on the measured values of a cp (from the one parameter fit) and c cp (from 
the two parameter fit). 



different numbers of charm hadrons may exhibit larger CP-violating effects of different signs which 
largely cancel out in the inclusive decay. In this case, the expected asymmetry takes the form 

A{t') = c cp sin(Am d t') - a cp sin 2 ( \ (16) 

where a cp is still related to Re(ee) as in equation [14, and c cp is an additional CP-violating parameter 
[^]. A second fit to the form shown in equation |l£, allowing the values of both a cp and c cp to vary, 
gave the results 

a cp = 0.002 ±0.055 
c cp = 0.026 ±0.027 

where again the errors are statistical only, and the correlation coefficient between the two parameters 
is 0.72. No evidence is seen for a significant sin(Am^t') term, and the value of a cp shifts by only 
—0.003 (5 % of the statistical error) with respect to the one parameter fit, showing that the effects of 
efficiency differences are not significant. 

The systematic errors on the values of a cp and c cp (from the one and two parameter fits respectively) 
are all small compared to the statistical errors. They are shown in Table |3|, and discussed in more 
detail below. 

Physics input parameters: The values of the B° lifetime r B o and oscillation frequency Amj were 
taken from Q and varied within the quoted errors. 

B° fraction: The fraction of B° events in the data sample depends on the production fractions of B s 
and Ab in bb events M and on the fraction of non-bb background, determined to be 13.7 ± 3.0 % 
from the data. 

Flavour tagging: The offsets in the jet and vertex charges used for tagging the B° production flavour 
were measured directly in the data, as described in Section |5|. The uncertainties in these offsets 
contribute to the systematic errors on a cp and c cp . 

If a time dependent CP-violating effect were present, it could change the value of the offset 
measured for the T-tag vertex charge Q v tx, as the vertex tagged sample (being biased towards 
long B° decay times) would contain an unequal mixture of B° and B° which do not have equal 
vertex charge offsets. With the offsets measured in Monte Carlo and a CP-violating effect of 
flcp = 0.05, this effect would shift the Q v t x offset by 0.0006 of the RMS width of the Q v t x 
distribution. This shift is much smaller than the statistical error on the vertex charge offset in 
the data, and does not present an important additional source of systematic error. 
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The production flavour mis-tag probabilities as a function of Qt and Qm were also tested in 
the data, and found to be correct to precisions of 2.6% and 10%, respectively (see Section ||). 
These uncertainties were translated into errors on a cp and c cp by scaling the values of Qt and 
Qm by ±2.6% and ±10% as in \27\ and repeating the fits. 



Reconstruction asymmetry: The fit assumes that the secondary vertex reconstruction efficiencies 
in the M-tag hemisphere are equal for B° and B° mesons. The track reconstruction asymmetries 
mentioned in section |5| may potentially introduce an efficiency asymmetry, since the vertex 
reconstruction requires a secondary vertex with at least three tracks, and the sign of the highest 
momentum tracks will be different for B° and B° mesons. However, the track reconstruction 
asymmetries in the Monte Carlo, which are somewhat larger than those in the data, lead to no 
significant efficiency asymmetry for reconstructing secondary vertices. 

b fragmentation: The proper time reconstruction depends slightly on the average energy of the b 
hadrons. This effect was assessed by reweighting Monte Carlo events and repeating the fit, as 



discussed in Section 6.4 



Time resolution: The fit is rather insensitive to the reconstructed proper time resolution, since the 
effects of CP violation are characterised by a time scale t ~ Tr/Am^ which is much larger than 
the average proper time resolution. The sensitivity was assessed by varying the width of each 
Gaussian in the resolution function by ±10%, varying the fraction of the wider Gaussian by 
±50% and by using an alternative resolution function parameterisation derived from a Monte 



Carlo sample with 10% degraded tracking resolution, as in Section 6.4. 



The fit was tested on Monte Carlo by introducing non-zero values of a cp and c cp , and checking 
that the fit correctly reproduced the input asymmetries. The fit errors were also verified to be correct 
by splitting the Monte Carlo input into several sub-samples and studying the distributions of fitted 
outputs. 

Including all systematic errors, the value of a cp was determined from the one parameter fit to 
be o cp = 0.005 ± 0.055 ± 0.013. The value of c cp was determined from the two parameter fit to be 
c cp = 0.026 ± 0.027 ± 0.015. The result for a cp can be translated into a measurement of Re(ee) using 



equation |L4J, and gives 

Re(e B ) = 0.001 ± 0.014 ± 0.003 
where the first error is statistical and the second systematic. 



7.3 Fit and results for (Ar/r) b 

The same data sample was used to measure the fractional difference between b and b hadron lifetimes. 
This was done by dividing the data into 20 bins of b/b hadron purity using the tagging variable Qi 
and simultaneously fitting all the reconstructed proper time distributions. The expected proper time 
distribution Fj(t) in bin j of the tagging variable Q2 (—1 < Q2 < 1) is given by: 

roc 

FM) = / PJt')R(t - t')dt' 
Jo 

where P(t') describes the true proper time distribution of the events and R(t — t') is the time resolution 
function. The true proper time distribution is given by: 

W) = (l-/ bg ) ( (1 + ( 9 92)j) -e-^b + (l-(^) l e -^ b } +/bg J W) + (1 _ /,) J- e -W 

where /b g is the fraction of non-bb background in the data sample and {Q2)] is the average value of 
Q2 in the tagging bin j. The parameters Tb and are the lifetimes of b and b hadrons, related to the 
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Source 


A((Ar/r) b ) 


Background fraction 


0.000 


Background lifetime 


0.000 


Flavour tagging offsets 


0.001 


Flavour tagging mis-tag 


0.008 


b fragmentation 


0.000 


Time resolution 


0.001 


Total 


0.008 



Table 4: Systematic errors on the measured value of (At/t)^. 



average lifetime T av and fractional difference (At/t^ by equation [L5|. The non-bb background was 
modelled by two components — a fraction f$ at t' = and the remainder having a lifetime r^ g , both 
distributions being convolved with the same resolution function as the signal. 

A two parameter x 2 fit was performed to the 20 reconstructed proper time distributions in 0.5 ps 
bins between 1 and 15 ps, fitting the values of (Ar/r)b and r av . The region below lps was excluded 
from the fit because the function Fj(t) does not give a good representation of the data in this region. 
This is because of the effects of the primary vertex on the resolution function R(t — t') discussed in 
Section 7.1. Such effects have to be taken into account when fitting the lifetime itself |HJ, but are not 
important when searching for a difference in b and b hadron lifetimes. 



In the fit, the background fraction was set to /b g = 13.4 ±3.0 % as in Section 7.2. The background 
parameters were taken from the Monte Carlo and set to fs = 0.87 and Tb g = 2.9 ps. 
The results of the fit were 

(^-pj = -0.001 ±0.012 

and r av = 1.500 ± 0.003 ps, where the errors are statistical only. The result for r av has large systematic 
errors of the order of 0.1 ps associated with the resolution function and should not be interpreted as 
a measurement of the average b hadron lifetime. 

The distributions of proper time in four different ranges of Q2, together with the fit results and 
expectation for (Ar/r)b = 0.2, are shown in Figure ||. The fit describes the data well, and no evidence 
for a difference between b and b hadron lifetimes is seen. 

The systematic errors on (Ar/r)b are small, as most uncertainties affect both lifetimes in the same 
way. They are summarised in Table |j. Most of them are similar to those of the CP-violation analysis 
described in Section |7.2j , and were evaluated in the same way. The uncertainty from the background 
lifetime was evaluated by varying the fraction of background at t' = from 82% to 92%, and by 
varying the lifetime of the other component from 2.4 ps to 3.4 ps. No significant effect was seen. 

The total systematic error on (Ar/r)b is ±0.008, and is dominated by the uncertainty on the 
flavour mis-tag rates. Additionally, it was checked that the fit was stable with respect to variations 
of the minimum and maximum time cuts. The fit was tested on Monte Carlo by reweighting so as 
to introduce variations between the b and b hadron lifetimes, and checking that the fit recovered the 
correct values of (Ar/r)b. The errors were checked by splitting the Monte Carlo into subsamples, as 
before. 
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Figure 9: Distributions of reconstructed proper time t (1 < t < 15 ps) in four bins of tag variable Qi- 
The data are shown by the points with error bars, and the prediction of the fit by the solid lines. The 
expected distributions for a 20% difference between b and b hadron lifetimes ((At/t% = 0.2) are 
shown by the dotted lines. 
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8 Summary and conclusions 

The lifetimes of the B + and B° mesons, and their ratio, have been measured using a technique based 
on reconstructed secondary vertices. From a sample of 2 390 221 hadronic Z° decays collected by the 
OPAL detector between 1993 and 1995, the results 

r B+ = 1.643 ± 0.037 ± 0.025 ps 
r B o = 1.523 ± 0.057 ± 0.053 ps 
r B+ /r B o = 1.079 ± 0.064 ± 0.041 

were obtained, where in each case the first error is statistical and the second systematic. These results 
are in agreement with other measurements from LEP, SLD and CDF |3|-[6| , and with the world average 
values of t b + = 1.65 ± 0.04 ps and r B o = 1.56 ± 0.04 ps [jij. The result for the B + lifetime is the most 
precise determination to date. 

Using a similar technique, an inclusive sample of b hadron decays has been used to search for CP 
and CPT violation effects. No such effects have been seen. From the time dependent asymmetry of 
inclusive B° decays, the CP violation parameter has been determined to be 

Re(e B ) = 0.001 ± 0.014 ± 0.003 . 

This result agrees with the OPAL measurement using semileptonic b decays: Re(e B ) = 0.002 ± 0.007 ± 
0.003 jlOfl , and is also in agreement with other less precise results from CLEO and CDF ||. 
The fractional difference in the lifetimes of b and b hadrons has also been measured to be 

/ At\ r(b hadron) — r(b hadron) 

= — i-r = -0.001 ±0.012 ±0.008 . 

V r J h r(average) 

This is the first published analysis to test the equality of the b and b hadron lifetimes. 
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